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Abstract 

A  real-time  thermoregulatory  model  was  developed  for  predicting  real-time  physiological  responses  of  workers  engaged  in  various  tasks  for 
prolonged  time.  The  unique  feature  of  the  present  model  is  primarily  on  metabolic  activity  inputs  derived  from  minimum  non-invasive  measures 
(i.e.,  heart  rate  and  ambient  temperature).  In  addition,  it  utilizes  individual  anthropological  characteristics  (height,  weight,  and  clothing)  as  an 
input  to  estimate  core  temperatures  (Tc).  The  model  was  validated  using  data  from  five  laboratory  studies  ( n  =  63)  with  varied  environments, 
clothing,  and  heat  acclimation  status.  Overall,  Tc  predictions  using  this  simplified  model,  corresponded  well  with  measured  values  (root  mean 
square  deviation:  0.05-0.31  °C). 

Published  by  Elsevier  Ltd. 

Keywords:  Real-time  modeling;  Heat  stress;  Heart  rate;  Core  temperature;  Air  temperature 


1.  Introduction 

Thermoregulatory  mathematical  models  have  been  increas¬ 
ingly  utilized  to  assess  thermal  strains  of  workers  without  risk, 
cost,  and  time  related  human  experiments.  Although  statis¬ 
tical  models  limit  simulation  conditions  based  on  model  as¬ 
sumptions,  mathematical  models  consisting  of  series  of  thermal 
physiological  equations  can  provide  simulations  by  generating 
various  combinations  of  environmental  and  operational  heat 
stress  situations.  Previous  thermal  regulatory  models  require 
body  core  temperature  (Tc)  as  an  input  [1-3],  which  is  an  im¬ 
portant  physiological  parameter  for  the  assessment  of  thermo 
physiological  strain  as  it  reliably  indicates  impending  injury 
[4,5].  However,  Tc  is  an  invasive  measurement  and  often  im¬ 
practical  and  undesirable  measurement  for  monitoring  during 
varied  activities  over  prolonged  periods  in  industrial  and  mili¬ 
tary  situations.  In  addition,  utilizing  previous  models  to  predict 
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physiological  responses  can  be  inconvenient  and  inaccurate 
for  workers,  especially  soldiers,  engaged  in  various  strenuous 
field  activities  for  prolonged  time.  For  instance,  some  models 
can  predict  physiological  responses  based  on  only  consistent 
environmental  and  operational  conditions  during  work  con¬ 
trolled  duties  [2,3].  Furthermore,  many  inputs,  some  of  which 
may  be  invasive,  are  required  for  previous  models  and  tend  to 
provide  more  inappropriate  or  missing  data  in  field  working 
environment  because  strenuous  activities  in  a  field  are  likely  to 
cause  sensor  malfunctions  [1,6].  Therefore,  this  new  suggested 
model  based  on  a  minimum  number  of  non-invasive  inputs 
was  needed  for  predicting  real-time  physiological  responses  of 
soldiers  during  various  prolonged  and  heat  stress  related  oper¬ 
ations.  The  purpose  of  this  paper  is  to  develop  the  new  model 
and  test  the  accuracy  of  the  model  predictions  with  a  focus  on 
Tc  using  available  data.  Although  the  model  in  this  study  was 
primarily  designed  for  military  personnel  who  are  usually  heat 
acclimatized  during  training  with  normal  or  low  levels  of  body 
fat  relative  to  the  general  population,  the  concept  of  the  model 
can  be  applicable  to  industrial  and  emergency  operations. 
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Comparisons  of  physiological  data  with  model  predictions 
provide  the  analytical  basis  necessary  to  characterize  model 
performance  and,  may  if  necessary,  identify  areas  of  further 
improvements. 

2.  Methods 

The  new  model  represents  the  human  as  two  active  physio¬ 
logical  compartments  (core  and  skin)  surrounded  by  a  passive 
clothing  compartment.  The  initial  temperature  in  each  compart¬ 
ment  was  set  lower  (ri0.3DC)  for  heat  acclimated  individuals 
[6].  Within  a  compartment  the  thermophysiological  properties 
are  assumed  to  be  uniform;  that  is  for  example,  Tc  has  the  same 
value  at  time  t  throughout  the  core  compartment  and  likewise 
for  skin  temperature.  All  metabolic  heat  production  ( M )  occurs 
in  the  core.  Some  of  the  core  heat  is  lost  directly  to  the  envi¬ 
ronment  by  respiration  with  the  remaining  flowing  to  the  skin 
by  conduction  and  skin  blood  flow.  Overall,  the  regulation  of 
Tc  in  the  model  is  achieved  primarily  by  controlling  blood  flow 
to  the  skin  and  sweating  rate  of  the  skin.  These  are  reasonable 
representations  and  simplifications  because  of  increased  blood 
circulation  in  warm  and  hot  conditions  [6].  The  basic  thermal 
regulatory  functions  including  heat  production,  transfer,  bal¬ 
ance,  and  thermal  strain  index  utilized  in  the  model  were  de¬ 
scribed  below. 

2.1.  Metabolic  rate  (M) 

A  unique  feature  of  this  model  is  the  use  of  real-time  esti¬ 
mates  of  M  (Watts)  derived  from  the  measured  heart  rate  (HR) 
and  environment  temperature  [7].  This  equation  was  utilized 
because  the  parameters  are  minimum  non-invasive  measures 
and  overall  good  correlation  between  O2  uptakes  and  HR  were 
demonstrated  [7,8 j.  The  M  equation  using  HR  and  environmen¬ 
tal  temperature  was  described  as  follows: 

M  =  [0.68  +  4.69(HRratio  -  1) 

-  0.052(HRratio  -  l)(ra  -  20)]58.1Ao  (1) 

where  HRratio — observed  HR  given  at  the  time/resting  HR 
of  the  individual,  and  Ta — ambient  temperature  in  °C.  Ad  is 
the  body  surface  area  (nr).  Cardiac  output,  indicated  by  HR, 
supplies  blood  and  oxygen  for  metabolism  and  also  blood 
flow  to  the  skin  for  thermoregulation.  For  a  given  metabolic 
activity,  the  HR  will  increase  with  increasing  environmental 
temperature  because  of  the  thermoregulatory  need  for  in¬ 
creased  skin  blood  flow.  The  prediction  of  M  shown  in  Eq.  (1) 
was  derived  from  HR  and  oxygen  consumption  measured 
over  a  range  of  metabolic  acdvity  and  environmental  temper¬ 
ature  [7],  The  equation  is  tested  for  the  conditions  includ¬ 
ing  20°C  ^  7'a  40°C,  1.2  <  HRratio  ^  2.1,  wind  speed 

%  1.25  ms-1,  dew  point  temperature  ^  20  °C. 

2.2.  Heat  transfer 

Heat  exchange  from  the  skin  to  the  environment  by  radiation 
and  convection  is  classified  as  dry  heat  exchange.  The  rate  of 


dry  heat  loss  (2dry,  Wm  2)  was  determined  by 

Gdry(=7sk  ~  7o)/(Rcit)  (2) 

where  Rcit  is  the  total  dry  thermal  resistance  between  skin  and 
the  environment  and  Ta  is  the  operative  temperature  of  the 
environment  [6],  T^  is  the  skin  temperature. 

Rclt  values  for  military  clothing  used  in  validation  process 
for  the  present  model  are  from  thermal  manikin  measurements. 
T0  is  the  average  of  ambient  air  (Ta)  and  mean  radiant  temper¬ 
atures  weighted  by  their  respective  heat  transfer  coefficients  as 
described  by  Kraning  and  Gonzalez  [6], 

The  evaporative  heat  loss  from  the  skin  (Qevap.  Wm-2)  is 
from  water  diffusion  through  the  dry  skin  and  from  evaporation 
of  areas  covered  with  sweat.  It  is  quantified  as  [9j: 

2evap  ==  (I  -  u00.06i;max  +  wE 

max  (3) 

where  the  maximum  rate  of  evaporative  heat  loss  from  the  skin 
surface  (Emax)  completely  covered  by  sweat  is  calculated  from 
saturated  vapor  pressure  (Torr)  of  water  at  the  skin  tempera¬ 
ture,  the  ambient  water  vapor  pressure  (Torr)  and  the  vapor  re¬ 
sistance  of  the  clothing  system  from  skin  to  the  surrounding 
environment.  Skin  wettedness  (w)  or  the  fraction  of  the  skin 
surface  covered  by  sweat  was  determined  from  the  rate  of  sweat 
secretion  to  the  rate  of  maximum  evaporation  from  completely 
wet  skin  [2]. 

The  heat  exchanged  via  respiration  (Qres.  W  m~2)  was  esti¬ 
mated  [9]  as 

Qtes  =  M/Afl(0. 0014(34  -  ra)  +  0.0023(44  -  Pa))  (4) 

where  PA  is  the  ambient  water  vapor  pressure  (Torr)  and  Ta  is 
the  ambient  temperature  (°C). 

2.3.  Heat  balance 

A  heat  balance  analysis  of  the  core  compartment  yields, 

M/Ad  —  2res  +  Qk  +  <2skbf  +  (Wc/A/,)cbt(dTc/dt)  (5) 

where  Ml  Ad  is  in  W  m-2,  Qles  is  respiratory  heat  loss,  Qk  is  pas¬ 
sive  heat  conduction  from  core  to  the  skin,  and  g^bf  represents 
heat  transported  by  the  blood  from  core  to  skin  compartments. 
The  last  term  (dTc/dr)  on  the  right  in  Eq.  (5)  represents  the  rate 
of  heat  storage  of  the  core  compartment  where  Wc  is  the  mass 
(kg)  of  the  core  ( ~  0.95  of  total  body  mass  in  this  case)  and 
specific  heat  of  body  tissue  (cbt)  is  0.97  Whkg-1  °C-1  [2].  The 
passive  heat  conduction  from  core  to  skin  (Qk)  was  estimated 
based  on  the  temperature  differences  between  skin  and  core  and 
the  constant  conductance  ( k )  value  of  the  tissue  between  two 
compartments  (5.28Wm-2oC-1)  [9],  Similarly,  the  heat  trans¬ 
ported  by  blood  flow  (Sskbf)  to  the  skin  is  determined  based 
on  skin  blood  flow,  which  is  modeled  proportional  to  changes 
in  core  and  skin  temperatures  from  their  set  point  temperatures 
(Tcset  =  36.8  °C;  I^set  =  33.7  °C)  and  specific  heat  of  blood  as 
a  constant  value  of  1 . 1 63  W  h  L- 1  °C- 1  [2] .  The  rate  of  Tc  (°C) 
change  ( dTJdt )  can  be  rearranged  by  Eq.  (5)  to  find  the  next 
Tc.  The  energy  balance  of  the  skin  compartment  (Wm-2)  is 
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Table  1 

The  summary  of  three  heat  study  conditions  used  in  this  study . 


Study# 

HSl 

HS2 

HS3 

HS4 

HS5 

Ta  (°C) 

27 

49 

30 

35 

40 

RH  (%) 

75 

18 

50 

45 

40 

Acclimation  status 

No 

Yes 

Yes 

Yes 

Yes 

Clothing 

HWBDU 

T-shirt  &  shorts 

HWBDU 

CPG 

T-shirt  &  shorts 

Activity  (W) 

277-350 

420-620 

510-840 

-450 

-410 

Duration  (min) 

—  170 

<  100 

110 

<40 

125 

N 

9 

5 

11 

8 

30 

Reference  (#) 

[111 

[12] 

[121 

1 131 

[141 

HWBDU — hot  weather  battle  dress  uniform  and  CPG — chemical  protective  garment. 


modeled  similarly  to  Tc  substituting  the  mass  (kg)  of  the  skin 
compartment  and  T,±,  instead  [2]. 

2.4.  Physiological  strain  index  (PSI) 

The  thermal  strain  index  called  PSI  was  utilized  based  on 
measured  HR  and  predicted  Tc  by  the  model  for  evaluating  on¬ 
line  heat  strain  levels  of  individuals  calculated  as  follows  [10J: 

PSI  =  5 (Tct  -  Tc 0)(39.5  -  Tc)~l 

+  5(HR,  -  HRoX  1 80  -  HR0)“ 1  (6) 

where  Tcq  and  HRq  are  the  initial  7'c  and  HR,  respectively.  Tct 
and  HR,  are  taken  at  a  given  time  t  during  the  heat  exposure 
[TO].  Unlike  other  thermal  indices,  PSI  is  simple,  easier  to  use, 
and  feasible  for  different  operational  and  environmental  con¬ 
ditions  [10].  The  PSI  consists  of  0-10  scale,  classifying  them 
into  five  thermal  categorical  states  (no  to  little,  low,  moderate, 
high,  very  high)  [10]. 

2.5.  Validation  data 

The  validation  of  the  model  was  conducted  using  data  from 
five  different  laboratory  studies.  In  all  cases,  participants  pro¬ 
vided  their  informed  consent  prior  to  participation  and  were 
free  to  withdraw  from  the  study  at  any  time. 

Heat  Study  1  (HS1):  Nine  unacclimated  volunteers  (8  males, 
1  female;  age:  23  ±4  [SD]yr;  height  174.2  ±  5.2  cm;  weight: 
73.4  ±  6.5kg),  wearing  the  US  Army  hot  weather  battle  dress 
uniform  (HVVBDU)  did  intermittent  exercise  in  a  warm-humid 
environment  (air  temperature  (Ta):  27  °C,  relative  humidity 
(RH):  75%RH)  for  170min  [11].  The  exercise  routine  con¬ 
sisted  of  walking  on  a  level  treadmill  at  1.34  ms-1  for  30  min 
followed  by  lOmin  of  rest.  Grand  mean  of  subjects’  HR  ranged 
between  74  and  106bpm  and  that  of  Tc  ranged  between  37.0 
and  37.5  °C  during  the  experiment.  Rectal  temperature  was 
measured  each  minute  [1 1  j. 

Heat  Study  2  (HS2):  Five  heat  acclimated  male  soldiers  (age: 
22 ± Syr,  height  176 ± 4cm;  weight:  71.1  ±  9.3kg)  wearing 
t-shirt  and  shorts  continuously  walked  at  420-620  W  in  a  hot 
dry  (49  °C,  18%RH)  environment  [12].  The  maximum  exer¬ 
cise  duration  was  100  min  and  subjects  continued  a  test  session 
until  they  voluntary  withdrew.  Grand  mean  of  subjects’  HR 


increased  from  83  to  157bpm  and  that  of  Tc  increased  from 
36.8  to  38.6 °C  during  the  experiment.  Tc  was  measured  with 
an  ingestible  telemetry  temperature  pill  [12]. 

Heat  Study  3  (HS3):  Eleven  volunteers  (10  males,  1  female; 
age:  21±8yr;  height:  175±7cm;  weight:  78.0±  11.4kg), 
wearing  HWBDU,  did  intermittent  exercise  in  a  warm-humid 
environment  (30°C,  50%RH)  for  llOmin  [12],  Subjects  are 
partially  heat  acclimated  by  natural  summer  acclimation  but  no 
official  acclimation.  The  exercise  routine  consisted  of  walking 
on  a  treadmill  (1.56ms-1, 1%  grade)  for  50min,  followed  by 
10  min  rest.  Grand  mean  of  subjects’  HR  ranged  between  82 
and  152bpm  and  that  of  Tc  ranged  between  36.8  and  38.2°C 
during  the  experiment.  Rectal  temperature,  measured  with 
an  ingestible  telemetry  temperature  pill  was  recorded  every 
10  min  [12]. 

Heat  Study  4  (HS4):  Eight  heat  acclimated  men  (age: 
23±6yr;  height:  176 ± 6cm;  weight:  76.0 ±  15.4kg),  wear¬ 
ing  chemical  protective  garments,  walked  on  a  treadmill  with  a 
4-9%  grade  at  1.56-1.65  ms-1  that  resulted  in  metabolic  effort 
about  450  W  and  oxygen  uptake  around  55%  of  V02max  until 
they  voluntary  withdrew  from  the  study  or  Tc  >  39.5 °C  [13J. 
Subjects  were  able  to  continue  exercising  for  <  40 min  under 
environmental  conditions  of  35°C/45%RH.  Tc  was  measured 
by  a  rectal  probe  [13]. 

Heat  Study  5  (HS5):  Thirty  heat  acclimated  male  volun¬ 
teers  (age:  19  ±  1  yr;  height:  180  ±  1  cm;  weight:  76.3  ±  2.8kg) 
wearing  t-shirt  and  shorts  walked  on  a  treadmill  at  ~4 1 0  W  in 
the  environmental  chamber  (40°C,  40%RH)  for  125  min  [14]. 
Grand  mean  of  subjects’  HR  increased  from  77  to  121  bpm  and 
that  of  7’c  increased  from  37.2  and  37.9 °C  during  the  experi¬ 
ment.  Rectal  temperature  was  measured  every  minute  utilizing 
a  thermistor  [14], 

Table  1  is  the  summary  of  heat  study  conditions  described 
above.  Depending  upon  the  individual  study,  HR  and  Tc  were 
collected  at  different  time  intervals.  Previous  studies  have 


suuvvu  uiv  iwiaunu) 

measured  Tc  [15,16]. 


2.6.  Statistical  analysis 

The  predicted  and  observed  Tc  for  each  individual  partici¬ 
pant  were  compared  using  root  mean  square  deviation  (RMSD) 
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method  as  described  in  Haslam  and  Parsons  [17].  The  RMSD 
was  used  to  quantify  the  average  difference  between  predicted 
and  observed  measurements  across  time  [17], 

3.  Results 

3.1.  Model  structure 

The  concept  of  the  model  described  above  was  integrated 
and  the  basic  operational  structure  of  the  new  model  was  sum¬ 
marized  in  Fig.  1.  The  model  can  use  individual  values,  group 
means  or  default  population  values  as  the  input  for  anthropo¬ 
logical  characteristics  (height,  weight,  or  clothing).  Then,  real¬ 
time  input  of  measured  HR  and  local  weather  (Ta,  wind  speed, 
RH,  and  radiant  load)  are  used  to  make  predictions  and  esti¬ 
mates  of  physiological  parameters  including  Tc,  sweat  rate,  ac¬ 
cumulated  water  loss,  T^,  metabolism,  and  thermoregulatory 
strain  [18].  The  shaded  area  indicates  an  input  variable.  A  box 
with  dotted  lines  displays  multiple  functions  of  the  model  for 
predicting  outputs  (a  box  with  solid  lines). 

3.2.  Validation 

Figs.  2-6  summarize  the  comparisons  of  the  mean  measured 
Tc  to  the  corresponding  mean  model  predicted  Tc  using  each  HS 
data.  Overall,  the  predictions  of  Tc  agreed  well  with  measured 
values  and  the  RMSD  ranged  between  0.05  and  0.31  °C. 

HS1:  In  a  moderate  environment,  subjects’  Tc  were 
<  37.5  °C  and  the  model  predictions  were  within  one  SD  of  the 
mean  measured  values  (Fig.  2).  The  mean  of  one  standard  de¬ 
viation  (SD)  between  subjects  across  the  time  was  0.19°C  and 
the  prediction  errors  for  HS1  were  small  (RMSD  =  0.10  °C) 
(Fig.  2). 

HS2:  Under  hot  dry  environment,  mean  Tc  of  volunteers 
raised  ~38.5°C  with  more  variability  between  subjects  across 


Fig.  1.  Basic  operational  structure  of  the  new  model. 


Fig.  2.  The  summary  comparisons  between  mean  measured  and  predicted 
core  temperature  (Tq)  from  Heat  Study  1  (27°C,  75%;  hot  weather  battle 
dress  uniform;  n  =  9). 


Experimental  time  (minutes) 

Fig.  3.  The  summary  comparisons  between  mean  measured  and  predicted  core 
temperature  (rc)  from  Heat  Study  2  (49°C,  18%;  t-shirt  and  shorts;  n  =  5). 


Fig.  4.  The  summary  comparisons  between  mean  measured  and  predicted 
core  temperature  {Tc)  from  Heat  Study  3  (30 °C,  50%;  hot  weather  battle 
dress  uniform;  n—  11). 


the  time  (0.23  °C)  than  HS1  (Fig.  3).  The  RMSD  for  the  model 
prediction  errors  was  0.31  °C.  Initially,  the  model  predicted  a 
faster  rate  of  increase  in  Tc  than  measured  Tc  for  these  heat 
acclimated  participants  exercising  49  °C  /18%RH.  However,  the 
predictions  improved  after  50  min  and  were  very  good  deriving 
the  second  hour  of  moderate  exercise  (Fig.  3). 

HS3:  During  the  intermittent  exercise  in  a  warm-humid  con¬ 
dition,  subjects’  mean  Tc  increased  to  38.2  °C.  The  variabil¬ 
ity  (SD)  was  also  increased  from  0.26  to  0.38  °C,  toward  the 
end  of  the  exercise  (Fig.  4).  The  model  prediction  followed  the 
patterns  of  Tc  responses  to  intermittent  exercise,  within  1SD. 
RMSD  was  0.18  °C  (Fig.  4). 

HS4:  During  the  intermittent  exercise,  subjects’  measured 
mean  Tc  raised  ~38.4°C  (Fig.  5).  More  variability  in  mea¬ 
sured  Tc  between  subjects  was  observed  toward  the  end  of  the 
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Fig.  5.  The  summary  comparisons  between  mean  measured  and  predicted 
core  temperature  (71)  from  Heat  Study  4  (35°C,  45%;  protective  garment; 
n  <  8). 
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Fig.  6.  The  summary  comparisons  between  mean  measured  and  predicted  core 
temperature  (Tc)  from  Heat  Study  5  (40°C,  40%;  t-shirt  and  shorts;  n  —  30). 

exercise,  indicating  different  physiological  responses  to  heat 
(Fig.  5).  The  prediction  of  Tc  in  this  study  agreed  with  mea¬ 
sured  values  within  1SD,  and  RMSD  for  the  model  prediction 
errors  was  relatively  small  (0.12 °C)  (Fig.  5). 

HS5:  Under  the  operational  heat  stress  (40  °C,  40%RH), 
subjects’  mean  Tc  increased  to  37.8  °C  during  the  exercise 
(Fig.  6).  The  mean  of  SD  between  subjects  across  the  time  was 
0.35  °C.  Tc  predicted  by  the  model  agreed  well  with  the  mea¬ 
sured  values  and  the  model  prediction  errors  were  very  small 
(RMSD  =  0.05  °C)  (Fig.  6). 

4.  Discussion 

The  new  model  was  developed  to  assist  with  real-time  phys¬ 
iological  monitoring  of  individuals  engaged  in  prolonged  work 
during  training  or  military  operations.  The  initial  analysis  of  the 
model  described  in  this  paper  yielded  encouraging  results  for 
predicting  the  thermal  status  of  individuals  or  groups  exposed 
to  various  heat  stress.  The  model  predictions  corresponded 
to  the  patterns  of  Tc  responses  to  different  exercises  includ¬ 
ing  intermittent  or  constant  exercise.  Predictions  for  Tc  were 
within  an  acceptable  range  ( ±  1SD)  among  five  laboratory 
studies  varying  different  environmental  and  operational  heat 
stress,  acclimation  status,  and  clothing  except  for  a  certain  ex¬ 
ercise  period  in  HS2.  Initially,  in  HS2,  the  prediction  of  Tc  was 


conservative  in  hot  dry  environment,  that  is,  it  tended  to  in¬ 
crease  faster  than  measured  Tc\  however,  after  approximately 
sixty  minutes,  Tc  predictions  were  in  good  agreement  with  mea¬ 
sured  values.  This  is  because  of  higher  HR  ratio,  resulting  into 
higher  M  estimates  than  actual  rates,  in  the  beginning.  Despite 
comparing  non-compensable  heat  stress,  that  occurring  with 
protective  clothing  (HS4),  predictions  represented  measured  Tc 
fairly  well  (RMSD  =  0.12°C).  Normally,  soldiers  deploying  to 
regions  known  for  extreme  heat  will  be  acclimatized  during 
training  for  maximizing  their  performance  and  prevention  from 
heat-related  injuries.  As  the  model  based  on  a  minimum  num¬ 
ber  of  non-invasive  measures  was  primarily  developed  for  heat 
acclimated  individuals,  an  acclimatization  will  provide  more 
accurate  predictions  of  soldiers’  physiological  status  than  unac¬ 
climated  individuals,  particularly  under  strenuous  environmen¬ 
tal  conditions. 


The  model  introduced  in  this  study  has  an  advantage  over 
other  heat  prediction  models  in  that  estimates  of  M  from  HR 
and  environmental  temperature  are  obtained  using  non-invasive 
methods.  Because  HR  and  oxygen  consumptions  tend  to  be 
linearly  associated  [7,8],  HR  is  useful  to  estimate  metabolic 
rates.  Other  methods  such  as  oxygen  consumptions  and  doubly 
labeled  water  are  alternative  to  measure  metabolic  rates;  how¬ 
ever,  they  can  be  destructive,  impractical  or  costly,  in  prolonged 
hours  offbee-range  working  operations  [8J.  In  addition,  invasive 
measures  (e.g.,  Tc)  may  interfere  with  performance  and  real¬ 
time  monitoring  of  soldiers  or  workers  engaged  in  long  hours  of 
training  in  a  wide  range  of  hot  environment.  The  model  predicts 
good  agreement  with  measured  Tc ,  which  is  useful  for  screen¬ 
ing  physiological  and  health  status  of  workers  who  may  be  re¬ 
quired  to  perform  different  tasks  over  many  hours  or  even  days. 
This  model  can  provide  additional  assistance  for  military  oper¬ 
ations.  For  instance,  recent  reports  from  Iraq  and  Afghanistan 
indicate  that  faster  access  to  medical  care  and  rapid  identifica¬ 
tion  of  probable  evacuations  would  save  more  soldiers’  lives 
[19,20].  One  current  strategy  to  shorten  the  time  between  in¬ 
jury  and  treatment  is  to  position  a  small  surgical  care  unit  near 
the  battlefield  [19].  The  model  can  contribute  to  this  effort  by 
forecasting  both  near  real-time  and  the  future  probability  of  in¬ 
dividual  soldier’s  physiological  status  and  by  directing  medical 
attention  to  the  individual  soldier,  unit  medical  personnel,  or 
command  elements. 


There  are  several  programs  underway  in  the  Department  of 
Defense  to  establish  methodologies  and  systems  for  the  accu¬ 
rate  “real-time”  measurement  of  physiological  status.  Although 
obtaining  “true”  values  of  environmental,  physiological  and  op¬ 
erational  conditions  are  desirable  for  optimum  assessment  of 
soldier  status,  the  reality  of  military  operations  or  other  stren¬ 
uous  effort  under  extreme  conditions  requires  consideration  of 
alternative  methods.  For  instance,  the  loss  of  sensor  signal, 
calibration  errors,  and  interference  between  sensors,  are  com¬ 
mon  occurrences  resulting  in  the  loss  of  physiological  measure¬ 
ments  and  accurate  health  assessments.  If  model  predictions 
are  complicated  by  the  requirements  of  many  physiological, 
operational,  and  environmental  parameters,  the  accurate  predic¬ 
tions  are  limited.  Strenuous  activity  may  cause  sensor  malfunc¬ 
tion,  or  make  it  difficult  to  repair  the  device  during  training  or 
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operational  activities.  The  present  model  provides  a  “rugged” 
alternative  to  more  elaborate  sensor  suites  by  utilizing  mini¬ 
mum  inputs  of  non-invasive  measurements,  and  can  easily  be 
adjusted  after  sensor  loss  or  to  recover  data. 

This  paper  has  described  die  promising  performance  of  a  new 
prediction  model  using  laboratory  data.  However,  to  further  as¬ 
sess  model  performance,  to  build  user  confidence  in  model  pre¬ 
dictions,  and  to  identify  possible  improvements  in  application 
and  hardware,  it  is  important  to  test  the  model  during  field  train¬ 
ing  involving  operational  and  environmental  conditions  that  are 
close  to  operational  conditions.  Although  the  present  version  of 
our  model  was  created  for  a  specific  application  using  a  highly 
trained,  heat  acclimated  population,  the  model  is  readily  adapt¬ 
able  to  population  with  variations  in  fitness,  heat  acclimation 
status,  body  mass,  age,  or  other  physical  characteristics  to  ex¬ 
tend  the  model’s  application  to  more  diverse  groups  of  popu¬ 
lations  and  environments. 

5.  Summary 

A  thermal  regulatory  model  was  developed  for  predicting 
real-time  physiological  responses  of  workers  engaged  in  var¬ 
ious  tasks  for  prolonged  time.  The  model  can  use  individual 
values,  group  means  or  default  population  values  as  the  input 
for  anthropological  characteristics,  real-time  input  of  measured 
HR  and  local  weather  to  predict  real-time  physiological  pa¬ 
rameters  including  Tc,  Tsk,  M,  water  loss,  and  thermoregula¬ 
tory  strain.  The  model  was  validated  with  various  operational 
(e.g.,  clothing,  activities,  acclimation  status)  and  environment 
tal  (e.g.,  27-49  °C)  conditions  from  five  laboratory  heat  studies 
(n  —  63).  The  Tc  predictions  were  overall  in  good  agreement 
(RMSD:  0.05-0.3 1°C)  with  measured  core  temperatures. 

Conflict  of  interest  statement 

None. 

Disclaimer 

The  investigators  have  adhered  to  the  policies  for  protection 
of  human  subjects  as  prescribed  in  Army  Regulation  70-25, 
and  the  research  was  conducted  in  adherence  with  the  provi¬ 
sions  of  32  CFR  Part  219.  The  opinions  or  assertions  contained 
herein  are  the  private  views  of  the  author(s)  and  are  not  to  be 
construed  as  official  or  as  reflecting  the  views  of  the  Army  or 
the  Department  of  Defense. 

References 

[!]  K.K.  Kraning,  R.R.  Gonzalez,  A  mechanistic  computer  simulation  of 
human  work  in  heat  that  accounts  for  physical  and  physiological  effects 
of  clothing,  aerobic  fitness,  and  progressive  dehydration,  J.  Therm.  Biol. 
22  (1997)  331-342. 

[2]  A.R  Gagge,  A.R  Fobelets,  L.G.  Berglund,  A  standard  predictive  index 
of  human  response  to  the  thermal  environment,  ASHRAE  Trans.  92 
(1986)  709-731. 

[3]  K.B.  Pandolf,  L.A.  Stroschen,  L.L.  Drolet,  et  ah.  Prediction  modeling 
of  physiological  responses  and  human  performance  in  the  heat,  Comput. 
Biol.  Med.  16  (1986)  319-329. 


[4]  D.  Amos,  R.  Hansen,  W.M.  Lau,  et  ah,  Physiological  and  cognitive 
performance  of  Soldiers  conducting  routine  patrol  and  reconnaissance 
operations  in  the  tropics.  Mil.  Med.  165  (2000)  961-966. 

[5]  K.B.  Pandolf,  R.F.  Goldman,  Convergence  of  skin  and  rectal  temperatures 
as  a  criterion  for  heat  tolerance,  Aviat.  Space  Environ.  Med.  49  (1978) 
1095-1101. 

[6]  K.K.  Kraning,  R.R.  Gonzalez,  Physiological  consequences  of  intermittent 
exercise  during  compensable  and  uncompensable  heat  stress,  J.  Appl. 
Physiol.  71  (1991)  2138-2145. 

[7]  L.G.  Berglund,  Heart  rate  as  an  indicator  of  metabolic  rate  in  hot 
environments,  in:  Alliance  for  Engineering  in  Medicine  and  Biology, 
Proceedings  of  30th  Annual  Conference  on  Engineering  in  Medicine  and 
Biology,  19,  1977,  Los  Angeles,  CA,  Arlington,  VA,  November  5-9, 
1977,  pp.  274. 

[8]  W.D.  McArdle,  F.I.  Katch,  V.L.  Katch,  Human  energy  expenditure  during 
rest  and  physical  activity,  in:  W.D.  McArdle,  F.I.  Katch,  V.L.  Katch 
(Eds.),  Exercise  Physiology:  Energy,  Nutrition  and  Human  Performance, 
Lippicott  Williams  and  Wilkins,  Baltimore,  1996,  pp.  150-165. 

[9]  P.O.  Fanger,  Thermal  Comfort,  MacGraw-Hill,  New  York,  1972. 

[10]  D.S.  Moran,  A.  Shitzer,  K.B.  Pnadolf,  A  physiological  strain  index  to 
evaluate  heat  stress.  Am.  J.  Physiol.  275  (1998)  R129-R134. 

[11]  W.R.  Santee,  L.G.  Berglund,  A.  Cardello,  et  ah.  Physiological  and 
comfort  assessment  of  volunteers  wearing  Battle  Dress  Uniforms  (BDU) 
of  different  fabrics  during  intermittent  exercise.  Technical  Report 
T06-06,  US  Army  Research  Institute  of  Environmental  Medicine,  MA, 
2006. 

[12]  S.N.  Cheuvront,  S.L  Montain,  D.A.  Goodman,  et  ah.  Evaluation  of  the 
limits  to  accurate  sweat  loss  prediction  during  prolonged  exercise,  Eur. 
J.  Appl.  Physiol.  101  (2007)  215-224. 

[13]  W.A.  Latzka,  M.N.  Sawka,  S .J.  Montain,  et  ah.  Hyperhydration:  tolerance 
and  cardiovascular  effects  during  uncompensable  exercise-heat  stress,  J. 
Appl.  Physiol.  84  (1998)  1858-1864. 

[14]  D.S.  Moran,  Y.  Heled,  L.  Still,  et  ah.  Assessment  of  heat  tolerance 
for  post  exertional  heat  stroke  individuals,  Med.  Sci.  Monit.  10  (2004) 
CR252-CR257. 

[15]  M.A.  Kolka,  M.D.  Quigley,  L.A.  Blanchard,  et  ah.  Validation  of  a 
temperature  telemetry  system  during  moderate  and  strenuous  exercise, 
J.  Therm,  Biol.  18  (1993)  203-211, 

[16]  C.  O'Brien,  R.W,  Hoyt,  M.J.  Buller,  et  ah.  Telemetry  pill  measurement 
of  core  temperature  in  human  during  active  heating  and  cooling,  Med. 
Sci.  Sports  Exerc.  30  (1998)  468-472. 

[17]  P.A.  Haslam,  K.C.  Parsons,  An  evaluation  of  computer-based  models  that 
predict  human  responses  to  the  thermal  environment,  ASHRAE  Trans. 
94  (1998)  1342-1360. 

[18]  M.  Yokota,  L.G.  Berglund,  Initial  capability  decision  aid  (ICDA)  thermal 
and  its  validation.  Technical  Report  T06-03,  US  Army  Research  Institute 
of  Environmental  Medicine,  Massachusetts,  2006. 

[19]  T.R.  Bilski,  B.C.  Baker,  J.R.  Grove,  R.P.  Hinks,  Battlefield  casualties 
treated  at  Camp  Rhino,  Afghanistan:  lessons  learned,  J.  Trauma  54 
(2003)  814-822. 

[20]  A.  Gawande,  Casualties  of  war — military  care  for  the  wounded  from 
Iraq  and  Afghanistan,  N.  Engl.  J.  Med.  351  (2004)  2471-2475. 

Miyo  Yokota  received  her  B.A.  (1988)  from  Kyoto  Women’s  University  in 
Japan,  M.A.  from  Western  Michigan  University  (1991),  and  Ph.D.  in  Biolog¬ 
ical  Anthropology  from  University  of  Tennessee,  Knoxville,  in  1997.  Cur¬ 
rently,  she  is  a  Senior  National  Research  Council  Associate  in  Biophysics  and 
Biomedical  Modeling  Division,  US  Army  Research  Institute  of  Environmen¬ 
tal  Medicine,  Natick,  MA.  Her  research  interests  include  human  variability 
in  body  composition  and  thermal  regulatory  modeling. 

Larry  G.  Berglund  is  a  registered  Professional  Mechanical,  Electrical,  and 
Bioenvironmenta!  Engineer  by  education  (University  of  Minnesota  (B.S., 
M.S.),  University  of  New  Haven  (B.S.),  Kansas  State  University  (Ph.D.)). 
His  professional  career  includes  Trane  Company,  Michigan  Technological 
University,  Pierce  Foundation  Laboratory /Yale  University,  Tohoku  University 
(Japan),  and  Kimberly  Clark  Corporation  before  joining  US  Army  Research 
Institute  of  Environmental  Medicine  in  2000.  His  focus  is  on  applications  of 
bio-thermal  modeling. 


M.  Yokota  et  al.  /  Computers  in  Biology  and  Medicine  38  (2008)  1187 -1193 


1193 


Samuel  Cheuvront  received  degrees  in  Nutrition  from  Indiana  University  of 
Pennsylvania,  USA,  in  1994  and  1996,  and  a  Ph.D.  in  exercise  physiology 
from  Florida  State  University  in  2000.  From  2000  to  2003,  Dr.  Cheuvront 
completed  a  Post  Doc  in  Environmental  Physiology  at  the  US  Army  Re¬ 
search  Institute  of  Environmental  Medicine  (USARIEM)  where  he  remains 
currently  as  a  Principal  Investigator.  His  research  centers  around  the  study 
of  environmental  and  nutritional  influences  on  Soldier  performance. 

William  Santee  received  a  B.A.  degree  from  Gustavus  Adolphus  College 
in  1970,  an  MS  from  the  University  of  Georgia  in  1978,  and  a  Ph.D.  from 
Indiana  State  University  in  1985.  He  has  been  at  The  US  Army  Research 
Institute  of  Environmental  Medicine  since  1985,  where  he  is  a  Research 
Physical  Scientist  in  the  Biophysics  and  Biomedical  Modeling  Division.  His 
research  interests  focus  on  biometeorology  and  the  biophysics  of  human  heat 
exchange. 

William  Latzka  received  B.S.  degrees  in  Biology  and  Health  and  a  M.S.  in 
Exercise  Physiology  from  St.  Cloud  State  University,  Minnesota  ,  USA,  in 
1980  and  i982,  and  a  Sc.D,  in  Applied  Physiology  from  Boston  University 
in  1996.  From  1996  to  2001,  Dr.  Latzka  completed  a  Post  Doc  in  Envi¬ 
ronmental  Physiology  at  the  US  Army  Research  Institute  of  Environmental 
Medicine  (USARIEM)  where  he  is  currently  the  Chief  of  Research  Support 
Division/Executive  Officer.  His  research  interests  focused  on  improving  Sol¬ 
dier  performance  in  extreme  environments. 


Scott  Montain  received  a  B.S.  degree  from  the  University  of  Wisconsin- 
LaCrosse,  USA  in  1981  and  a  Ph.D..  from  The  University  of  Texas  at 
Austin,  USA,  in  1991.  Since  1991,  he  has  been  with  the  US  Army  Re¬ 
search  Institute  of  Environmental  Medicine,  Natick,  MA,  USA  where  he  is 
a  Research  Physiologist  in  the  Military  Nutrition  Division.  His  research  fo¬ 
cuses  on  nutritional  strategies  to  sustain  performance  during  environmental 
stress. 


Margaret  Kolka  received  her  B.S.  (1974)  and  M.A.  (1976)  degrees  from 
Michigan  State  University  and  her  Ph.D.  in  Exercise  Physiology  in  1980 
from  Indiana  University.  She  joined  Thermal  Physiology  Division,  US  Army 
Research  Institute  of  Environmental  Medicine  (USARIEM),  Natick,  MA  as 
a  Research  Physiologist  in  1983  and  was  Chief,  Biophysics  and  Biomedical 
Modeling  Division,  USARIEM,  from  2004-2005. 


Daniel  Moran  received  his  Ph.D.  (hon.)  from  the  Sadder  School  of  Medicine 
in  Tel  Aviv  University  (Israel)  in  1993.  Since  1990,  he  has  been  with  the 
Heller  Institute  of  Medical  Research,  where  he  acts  as  the  Director  from 
2000.  His  research  interests  focus  on  physiological  evaluation  of  human 
performance  during  various  environments  and  stresses,  non-invasive  methods 
for  physiological  monitoring,  and  biomarkers  of  bone-remodeling. 


